Topological Response Theory of Abelian Symmetry-Protected Topological Phases in 

Two Dimensions 



Meng Cheng^ and Zheng-Cheng Gu^'^ 

''Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park, MD 20742 
^Institute for Quantum Information, California Institute of Technology, Pasadena, CA 91125, USA 



o 

(N 
D 

O 
(N 



I 

o 
o 



> 
m 
o 

oo 

O 

cn 



Department of Physics, 



California Institute of Technology, Pasadena, CA 91125, USA 
(Dated: February 21, 2013) 



We systematically characterize the symmetry-protected topological (SPT) phases with finite 
Abelian symmetries that can be described by Chern-Simons field theory. We propose a topological 
response theory to uniquely identify the SPT orders, which allows us to obtain an (almost)complete 
characterization of bosonic SPT phases with any finite Abelian symmetry group. We then compare 
our results with the recently proposed group cohomolgy approach and discuss the limitation of the 
Chern-Simons theory approach. We show that even for finite Abelian symmetry, there exist bosonic 
SPT phases beyond the Chern-Simons theory framework. We also apply the theory to fermionic 
SPT with Zm symmetry and find the classification of SPT phases depends on the parity of m: for 
even m there are 2m classes, m out of which is intrinsically fermionic SPT phases and can not be 
realized in any bosonic system. Finally we propose a classification scheme of fermionic SPT phases 
for any finite, Abelian symmetry. 



PACS numbers: 

Introduction. The classification of gapped quantum 
phases [ij of matter has led to the conceptual distinc- 
tion between phases with long-range entanglement (LRE) 
and those with short-range entanglement (SRE). Differ- 
ent patterns of LRE describe quantum phases with differ- 
ent intrinsic topological (IT) orders, which are typically 
characterized by their topological ground state degener- 
acy on torus and fractionalized excitations. In contrast, 
SRE states have none of these features and are indis- 
tinguishable with the trivial product states (e.g. atomic 
insulators) from their bulk excitations. However, when 
symmetry is taken into account, finer structure emerges 
within the SRE states. A notable example is the recently 
discovered time-reversal-invariant topological insulators 
in two and in three dimensions 043 ■ These states are 
topologically nontrivial only if time-reversal symmetry 
is not broken (explicitly or spontaneously) and thus are 
called symmetry-protected topological phases @). Al- 
though SPT phases are seemingly featureless in the bulk, 
they do support gapless boundary excitations protected 
by the symmetry which clearly distinguishes them from 
trivial product states. 

The knowledge of SPT phases have been greatly ex- 
panded recently, including a complete classification of 
gapped quantum phases in one dimension and gen- 



eralizations to higher dimensions |13l- isj. In the gen- 
eral scheme, d-dimensional SPT phases in bosonic sys- 
tems with symmetry group G can be systematically con- 
structed from the {d + l)-th cohomology of G, denoted 
by H'^'^^{G, U(l)). Later on the construction is also gen- 
eralized to fermionic systems, using the so-called group 
super-cohomology approach [16|. These powerful con- 
structions are applicable to any symmetry group in arbi- 
trary spatial dimensions and have the advantage of giving 



directly the fixed-point ground state wavefunctions, as 
well as exact solvable parent lattice Hamiltonians. How- 
ever, the wavefunctions and parent Hamiltonians are gen- 
erally quite complicated and it is not easy to access the 
low-energy universal properties of these phases, e.g., the 
edge properties. 

Recently, it is proposed fl7-21| that a large family of 
SPT phases in two dimensions protected by Abelian sym- 
metries can be effectively described by Abelian Chern- 
Simons theory. The classification is then derived from the 
equivalence classes of K matrix and the (perturbative) 
stability of edge theory. However, a (non-perturbative) 
bulk argument is missing and the equivalent classes of 
K matrix is also quite difficult to compute for arbitrary 
Abelian symmetry group. On the other hand, the under- 
lying connection with (super)cohomology class is unclear 
as well. 

In this work we present a unified treatment of both 
bosonic and fermionic SPT phases in two dimensions, 
with finite, unitary Abelian symmetries. We generalize 
the proposal of Levin and Gu 22| , originally in the con- 
text of bosonic SPT phase protected by Z2 symmetry. 
The idea is to probe the bulk properties of SPT phases by 
coupling to a gauge field with gauge group G and examing 
the braiding statistics of gauge fluxes, as a unique "topo- 
logical response" of the SPT phases. This allows us to 
obtain an (almost) complete characterization of bosonic 
SPT phases with finite Abelian symmetry group that can 
be described in the Chern-Simons theory framework. We 
also study fermionic SPT phases with symmetry and 
interestingly find an even-odd effect: the classification is 
Z„i for odd m and for even m. Finally we obtain a 
(minimal) classification of fermionic SPT phases for any 
finite, Abelian symmetry. 
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Topological Response Theory of SPT Phases. First we 



review the Chern-Simons field theory of SPT phases [17, 
[isj . They are described by a multi-component Abelian 
Chern-Simons theory, which in its most general form is 
given by the following Lagrangian: 



CR = — e*^ a 



47r 



(1) 



Here K is a A'^ x iV integer symmetric matrix. Notice 
that for bosonic systems, all diagonal entries of K must 
be even while at least one of the diagonal entries is odd for 
fermionic systems. Since a SPT phase can not have any 
intrinsic topological order, we require detK = ±1 and 
there should be no chiral central charge on the boundary 
which means that the signature of K matrix must vanish. 

We will now assume that the system under consid- 
eration has a global, on-site symmetry group G. We 
only consider unitary, finite Abelian group, which gen- 
erally can be written as G = Z^i x x ■ • • Z^^, where 
TOfc > 1. We assume that the Z^^ subgroup is generated 
by gj with g^ ^ = 1. To understand the consequences 
of the symmetry, we couple the SPT phase to a gauge 
field taking value in the gauge group G. This is achieved 
by writing down minimal coupling for the matter fields 
in the SPT phases. The minimal coupling is completely 
specified by the irreducible representations carried by the 
matter fields. Since the symmetry is Abelian and finite, 
it amounts to assign Z^^ xZ^a x • • • Z^^. gauge charges to 
each matter fields. Therefore each matter field is associ- 
ated with gauge charges qa, valued in {0, 1, . . . , — 1} 
and the subscript a refers to different symmetry sub- 
groups. 

To facilitate the continuum field theory formulation, 
we view the group G as a subgroup of U(l)'"'. This per- 
spective allows us to introduce k external U(l) gauge 
fields A'^ ,a = 1 , . . . , fc that minimally couple to the mat- 
ter fields in the SPT phases and then by introducing k 
Higgs fields ipa with charge we can obtain a contin- 
uum version of discrete gauge theory as the symmetry 
broken phase of the U(l)'^ Higgs theory. The presence 
of Higgs condensation results in quantization of gauge 
fiuxes. 

After we gauge the global symmetry, the minimal cou- 
pling between the matter fields in the bulk and the ex- 
ternal gauge fields reads: 



-^^coupling — ^ ^ QqJ^ 



(2) 



The full theory is given by £ = £cs +'Ccoupiing- We then 
substitute j'" = ^e^^^^duB-x and integrate out internal 
gauge fields a, yielding an effective theory for the external 
fields A": 



We denote Ka/j = q^K~^q^ in the following discussion. 
We have therefore derived an Abelian Chern-Simons- 
Higgs theory as the effective "response theory" for the 
gauged SPT phases. Physically, as long as the symme- 
try G remains unbroken, we can always gauge the SPT 
phases. We propose that the braiding statistics of gauge 
fiuxes in the resulting theory fully characterize the SPT 
phases. 

Let us expose the general structure of the fluxes statis- 
tics starting from the simplest case G — Im- Then 
we just have a single Chern-Simons term at level K = 
q^K~^q. An "elementary" vortex in the Higgs fleld en- 
closes ^ gauge flux, as a result of flux quantization. Due 
to the screening of the Higgs fields, charge excitations 
in the Higgs field are all bosonic and gauge fluxes can 
form bound state with the Z^ charges, which are called 
"dyon" in the field theory context. The exchange statis- 
tics of the vortices can be readily read off from the Chern- 
Simons term ^AdA which binds Q = ir- charge to the 
vortex. For K not a multiple of 2to, this "fractional" 
charge can not be screened by the Higgs condensate. Due 
to Aharonov-Bohm effect, the flux has exchange statistics 
9 = The same result can be derived more carefully 
using charge-vortex duality. As we have argued, 9 (or 
equivalently, K ) serves as a unique "topological invari- 
ant" of the SPT phases with Z™ symmetry. However, we 
are still not done yet. A full classiflcation requires de- 
termining the equivalence classes of 9, which should be 
derived from the anyon spectrum in Chern-Simons- Higgs 
theory 23]. A heuristic, physical argument is the follow- 
ing: In deriving 9 we only consider a Higgs vortex. In 
principle there can be Z„i charges bound to it. Assum- 
ing q elementary charge is attached to the vortex, 
the composite object has statistical angle 9 + . This 
implies that the exchange statistics of flux as a "topo- 
logical invariant" for SPT phases is defined module 
or K is defined module 2m. Two SPT phases whose cor- 
responding statistical angles differ by an integer multiple 
of 1^ should be considered as being equivalent, since the 
resulting gauge theory is topologically indistinguishable 
(up to relabeling of dyons) . Thus the equivalance classes 
of K is at most Z2m- The Abelian group structure be- 
tween SPT phases is rather obvious since 9 is additive. 

However, this is not the end of the story. The statis- 
tics of the underlying particles, bosons or fermions, make 
a big difference. We first consider the simpler case of 
bosonic systems. Let us show that K for bosons must 
be even. This can be seen as follows: we write K — 
qi{K-')ijqj - 2E/<,7'Z/(K-i)/jg,7 + E/9l(K-i)//. 
Since all the diagonal entries of are even, it is then 
easy to see that K is even. Thus the classification is 
reduced to Z.^. They can be realized by the following 
{K,qJ: 
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£eff = — e^'^^A^q^K-iq^^.A^ + Cm^^,[^a,A^]. (3) 
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K = 



1 

1 2n 



,n = 0,l,.. 



1. 



(4) 
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— lira , we move to the next level 
X Z„ . In this case, we have to 



It corresponds to = ^ 
Having understood G 
of complexity G ~ Z„i 
consider an U(l) x U(l) external gauge fields and the 
response theory is fully characterized by a 2 x 2 matrix 
K. Repeating our argument above, the two diagonal ele- 
ments determine the exchange statistics of the two types 
of gauge fluxes corresponding to the two subgroups and 
give Z,„ X Z„ classiflcation. The off-diagonal element in- 
troduces a new ingredient, the braiding statistics between 
the two gauge fluxes, denoted by Q\2: B\i = ■ To 

determine the equivalence classes of 0\i^ we notice that 
we can attach charge to Z„ flux or vice versa and 
the mutual statistical angles are changed by integer mul- 
tiples of and respectively. Therefore we identify 
the equivalence relation as 



27rfci 27rfc2 



,ki,k2 & Z. 



In terms of K12, 



K12 = K12 + kill + k2m. 



(5) 



(6) 



The set of integers generated by kin + k2m, fci, ^2 G Z are 
simply the integer multiples of (m, n). Therefore we con- 
clude that that K12 is defined module (m, n), which gives 
additional „) classes. The complete classification is 
thus X Z„ X Z(„ „). 

What is the physical meaning of the mutual braiding 
statistics in the SPT phases? To give a concrete exam- 
ple, consider a SPT phase characterized by Ku — K22 = 
0, K12 = I. One can easily write down the following 
Chern-Simons theory: 



K 



1 

1 



,qi 



,q2 



(7) 



So the two "dual" degrees of freedoms in the SPT phases 
carry the two subgroup symmetries respectively. Other 
classes can be understood in a similar way. We however 
note by pass that 2 x 2 K matrix is not sufficient to 
realize all the SPT phases in this case if (m, n) > 1. 
For example, when m = n = 2, 7 out of the 8 SPT 
phases can be realized by 2 x 2 K matrix but the one 
with Kii = K22 = 2,i^i2 = 1 needs minimally 4 x 4 K 
matrix. In general, all bosonic SPT phases with Z,„ x Z„ 
symmetry are realized by 



K 



l2x2, Ql — 









1 







1 


,q2 = 


1 


w 







(8) 



It corresponds to K 



This is simply stacking 



'2p I 

up a layer of SPT phase defined in ([7]) with another layer 
which is essentially a Z„ (Zm) SPT phase. 



It is important to emphasize that the classification of 
SPT using the braiding and exchange statistics of fluxes 
is not equivalent to classifying the intrinsic topological or- 
der described by the corresponding dynamic gauge theory 
emerging at low energy. To classify the intrinsic topologi- 
cal order of the emergent gauge theory, only the braiding 
statistics is concerned and we do not keep track of the 
charge and flux carried by the quasiparticles. For exam- 
ple, when classifying by the intrinsic topological order we 
are free to exchange the labeling of gauge fluxes which 
is not allowed in the response theory. The difference is 
mostly manifested when m = n, since then the two gauge 
fluxes are "indistinguishable" , in the sense that they have 
the same spectrum of exchange statistics. For example, 
when TO = n = 2, the response theory gives Z^ clas- 
siflcation but the 8 theories collapse to only 4 different 
topological order (see the Supplementary Material [24 1 
for more details). In general, whenever {m,n) > 1 the 
two perspectives give different results. This issue has not 
been clarified in some related works 25, 261 in the con- 



text of classifications of symmetry-enriched topological 
orders, however, it is manifested in our approach. 

We are now well prepared to generalize the above pic- 
ture to arbitrary finite Abelian group G with k ^ 2 
generators. The K matrix of the response theory has 
k diagonal elements, yielding Z^^ x Z^^ x • • • Z^.^. clas- 
sification. The ^^^""'^^ off-diagonal elements determine 
the braiding statistics between the gauge fluxes and 
the classiflcation is Y[i<j'^(mi,mj)- Therefore, we con- 
clude that for a given symmetry group G = HiLi^mi, 
the Abelian Chern-Simons theory construction can give 



<i<j<k ^(mi,mj) 



classification of possible 



SPT phases. 

One may wonder whether this is the complete classifi- 
cation of SPT phases with symmetry group G. A totally 
independent approach, based on discrete topological 9- 
term and group cohomology taken by Chen et. al. shows 
that for a two-dimensional bosonic system with symme- 
try group G, the classification of bosonic SPT phases is 
given by the third cohomology of G with U(l) coefficients 
H^{G,U{1)). For G = IlLi^m.' one finds Q 

i7^(G,lLJ(l)) = JJ^Zm^ ]^Z(„j^^m^.) Y[ ^(m.,mj,mfe)- 
i i<j i<j<k 

(9) 

Clearly, when fc > 3 the Abelian Chern-Simons theory 
approach fails to capture the classes of SPT phases as- 
sociated with Z(„j. „j. „jj_-). In fact, it has been known 
that the corresponding group cocyles can not be realized 
as Abelian Chern-Simons theory. Physically, this is be- 
cause the gauge fluxes in the topological response theory 
carry non- Abelian statistics [23|. 

Fermionic SPT Phases. We now study fermionic SPT 
phases protected by a symmetry group G = "Lm- The 
conservation of total fermion parity is considered as a 
physical constraint rather than a symmetry. We will 
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show quite interestingly, the classification exhibits an 
even-odd effect: for odd m the classification is Z^, while 
for even to, it is 1j2m- The Z„ classes for odd m are 
nothing but the classes for bosonic SPT phases. 

The main difference between fermionic systems and 
bosonic ones is that the "identity" particle is a fermion 
which itself has a statistical angle tt. This implies that 
when we consider the equivalence classes of the braiding 
statistics of fluxes in the gauged theory, we are free to 
attach fermions to them and the statistical angle can be 
changed by n [27|, |28[ . We already see there is a differ- 
ence between even and odd m, since Z,„ gauge theory 
with even m contains a dyonic excitation with fermionic 
statistics while for odd m this is not the case. 

We first consider m odd. The Z,„ gauge theory itself 
allows one to change the exchange statistics of elementary 
gauge fluxes by = 0, 1, . . . , m — 1. By attaching 

fermions, the set of phases is fractionalized to = 
0, 1, . . . , m — 1. Given two SPT phases characterized by 
Ki and K2, the equivalence relation can be derived: 



ttKi 'KK2 



Trp 



(10) 



which implies that K = K + mp, in constrat io K = 
K + 2mp for bosonic SPT phases. We then easily see 
that there are only m different classes (We allow both 
even and odd K in the fermionic case). Therefore, all 
fermionic SPT phases with Z„i symmetry can be identi- 
fled with the bosonic ones. A physical interpretation is 
the following: bosonic SPT phases can be realized in the 
composite boson limit of any fermionic systems, in which 
pairs of fermions form bosonic molecules. This fact can 
also be derived directly (see the Supplementary Material 
for the derivation 2J]). 

Now we turn to the case when to is even. As we have 
elaborated, the Z^ gauge theory allows the statistical 
angle of fluxes to be changed by . Attaching fermions 
does not do much in this case since for p = to/2 we al- 
ready have a fermionic quasiparticle. The equivalence 
relation is exactly the same as the one for bosonic SPT 
phases, but now we allow both even and odd K. When 
K is odd, the corresponding SPT phases are intrinsically 
fermionic and can not be realized in any bosonic system. 
We are then led to the conclusion that the classification 
for to even is 1,2m ■ We notice that the to = 2 case receives 



a lot of attention recently |29l - l32| since it serves as a nice 



example of "collapse" of classification of non-interacting 
fermions (which is Z in the present case) when interac- 
tions are taken into account (Zg). The Chern-Simons 
field theory approach gives rise to a Z4 classification. 
The missing four classes are those with unpaired Ma- 
jorana edge modes which are beyond the Chern-Simons 
field theory description. 

Having worked out the classification for G = Z™, it is 
straightforward to generalize to arbitrary finite Abelian 
group G = Yia ^rn^ ■ Wc noticc that the underlying 



fermionic statistics does not effect mutual braiding prop- 
erties, so the additional classes due to nontrivial mutual 
braiding statistics between gauge fluxes in different con- 
jugate classes are stfll classified by JJ^^j 1{mi.mj)- Thus 
we obtained a classification of fermionic SPT based on 
Chern-Simons theory as Yli Im* Y[i<j ^{mi,mj) where to* 
is defined as to* ~ m for odd to and 2to for even to. 

Edge theory for fermionic SPT phase. Having estab- 
lished the classification of fermionic SPT phases with Z^ 
symmetry, we move on to study the edge theory and its 
stability in more detail. Let us consider the simplest one 
of Zm(for even to) fermionic SPT phases described by 
the following K matrix and q vector: 



K = 



1 

-1 



,q = 



(11) 



It is well known that the Chern-Simons theory im- 
plies existence of gapless edge states, whose effective La- 
grangian can be derived from gauge invariance princi- 
ple |33|: 



-Ccdgc = -^idt4'iKijdj:(t)j - 
47r 

with symmetry transformation: 



x(Pj) 



^ J, 27r 
^ H q 

TO 



(12) 



(13) 



To include interaction effect, it is convinient to switch to 
a non-chiral basis 4>i— ip — 9,(f>2 = (p + 9. We then have 
a generic Luttinger liquid model of the gapless edge: 



H 



dx ■ 



27r 



(14) 



Here u is the charge velocity and K the Luttinger pa- 
rameter. One can see that the Z^ transformation acts 
on the non-chiral bosonic fields as ip — > (p+^,d ^ + ^- 
To understand the stability of the edge theory, we add 
the leading perturbations cos 2mip and cos 2m9 allowed 

2 TV' 

by Zm symmetry. They have scaling dimensions ^ — 

2 

and ^respectively. So demanding that these two per- 
turbations are irrelevant, we find the stable region is 
<K 

2 

Conclusion and Discussion. In conclusion, we have 
systematically investigated the SPT phases with Abelian 
finite group symmetry in the Chern-Simons field theory 
framework. We develop a topological response theory to 
classify the SPT phases by gauging the symmetry group. 
By carefully studying the equivalence classes of the braid- 
ing statistics of gauge fiuxes, we characterize all possi- 
ble bosonic SPT phases that can be realized as Abelian 
Chern-Simons theories. We also compare our approach 
with the group cohomology theory approach and discuss 
the limitation of if-matrix approach. Indeed, the topo- 
logical response theory describes a non-perturbative ef- 
fect in the bulk of SPT phases and provides us a unique 



5 



way to identify different SPT phases. Finally, we extend 
the classification to fermionic SPT phases. For the sim- 
plest symmetry group G = Z„, we find the classification 
of fermionic SPT phases has an intriguing even-odd de- 
pendence on m. We then generalize the classification to 
arbitrary finite Abelian groups. We also discuss the edge 
stability of those intrinsic fermionic SPT phases. For 
future studies, it would be very interesting to develop a 
topological response theory to describe those SPT phases 
with non- Abelian flux statistics. On the other hand, the 
concept of topological response theory is also very use- 
ful for the classification of symmetry enriched topological 
order. 
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Supplementary material 



VORTEX DUALITY IN ABELIAN CHERN-SIMONS-HIGGS THEORY 



In this section we review the Abelian Chern-Simons-Higgs theory using charge-vortex duality transformation. The 
Lagrangian density of the theory is given by 



a 

The amplitudes of the Higgs fields are fixed by the potential energy. So we write (pa = 

Aff = ^e^-^A'-^K„pd,Af' - ^ Y(m„A^ - 8^9^^. 

a 

First we perform the Hubbard-Stratonovich transformation and write 

Aff = ^e^^'-'A-^K^ffd^A^ - ^ [^e^ - ^^{rn^A-, - d,e^) 



(SI) 



(S2) 



(S3) 
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Then decompose the phase field as 9a = rja + Ca where rja is the smooth part of the phase fluctuation and Cq is the 
singular (vortex) part. Integrate out the smooth part of the phase fields 77^, we obtain the constraint d^S^'^ = 0, which 
can be resolved as = 
The dual representation 

C = -^£^^^A"A-„^a,A^ - !^eM->A"a,6^ + . . . . (S4) 

47r ZTT 

We then integrate out the A fields to finish the duality transformation: 

/:duai = ^e^'-^b^mK-im^a.bA. (S5) 
From the dual action we can directly read off the self and mutual statistics of fluxes. 

If we treat the A fields as fully dynamical, the above action can be written as a doubled Chern-Simons theory: 



47r^ ' ' \M. 

Here M is the diagonal matrix diag(mi, 777,2, . . . ). We use the wedge product to simplify the expression. 

INTRINSIC TOPOLOGICAL ORDER AND CLASSIFICATION OF TOPOLOGICAL RESPONSE 

We elaborate on the nonequivalence between the classification of the gauged SPT phases as a "topological response" 
theory and the classification of the intrinsic topological order, using the example of the gauge group G = Z2 x Z2, 
which is also the symmetry group of the SPT phase. The emergent gauged theory is given by the following K matrices: 

02x2 22x2 \ Ik I 



The possible choices of fc,p and / yields all the H^{'L2 x Z2,U(1)) — Z2 cohomology classes. 

Here we emphasize that the classification of the response theory is not completely equivalent to the classification 
of the intrinsic topological order defined by the topological gauge theory (|S7|) where all the gauge fluxes are regarded 
as dynamical deconfined objects. In the latter case, two (Abelian) topological phases are equivalent as long as they 
have the same quasiparticle braiding matrices (so-called T and S matrices) , regardless of how the quasiparticles are 
labeled. This kind of equivalence relation is nothing but the §L(A^, Z) equivalence for the K matrix. In fact, the 8 
classes listed in (|S7I) reduce to only 4 under generic SL(4, Z) equivalence. In terms of the K matrices, there are only 
the following four different intrinsic topological orders: 

0\ (2 1 
2 

oyl^o 2) ^ \Q 2 

1\ (2 l\ fO 1 

1 oj I 1 0/ I 1 2 



(S8) 



Here ^ denotes the equivalence of the derived intrinsic topological orders. 

In the response theory, we are not allowed to permute the gauge fluxes from different subgroups of the symmetry 
group since they correspond to different physical symmetries. 

EQUIVALENCE BETWEEN FERMIONIC AND BOSONIC SPT PHASES WITH Z„ SYMMETRY 

We demonstrate directly that when 777, is odd all Z„j fermionic SPT phase are equivalent to bosonic ones. Let us 
consider the bulk Chern-Simons theory of Z™ fermionic SPT phase 



7 



Here g denotes the generator o£ the symmetry. We then add a trivial phase given by K = ax with trivial 
symmetry transformation on the edge bosons = where p ^Ij. We pack the whole system into a 4 x 4 K given 

by K = cr^ ® l2x2- 

We then perform the following §L(4, Z) transformation 



W 



/I 1 \ 
10 

-10 1 

VI 1 -ly 



(SIO) 



Under W the K matrix becomes K 



. So the first two components are describing fermionic systems and 



a,_ 
a^^ 

the last two bosonic ones. We denote the symmetry vector for the collectively as = {qi, q2,P,p)^ ■ The edge modes 
in the new basis are denoted by 



0: 



f4\ 

J 



W2I 



(p2 



(Sll) 



Under the §L(4, Z) transformation W, the vector q^ ^ W ^q^ = {qi — 2p, q2, 2p, qi — p)'^ . 

If qi — q2 is even, we let p = Si:^ and then q^ = ((72, 92, 9i — 92, ^^"5^)^. Thus the fermionic SPT is equivalent to 



a bosonic one with K 



(91 



q2,H^r. 



If qi — q2 is odd, it seems that p = being a half integer, is not physical. Here the fermionic nature plays a 



crucial. This is most easily understood from the edge modes. The edge modes (j) 
symmetry as 



h,4'2)'^ transforms under the 



Ug<PUl=<P- 



2np 
m 



(S12) 



When p is a half integer, we have 



U^4>{Ulr = <t> + 2T:p 



(p + n 



(S13) 



which is projectively the identity in a fermionic system. 

We now turn to the corresponding bosonic SPT. We notice that 





qi + q2 



(S14) 



which is not consistent with ?7™ = 1 in the bosonic case. Again the identity transformation in a fermionic system can 
be realized projectively: 



- TT 



(S15) 



which means that 



(S16) 



is also an identity transformation. So wc can freely add 



to (\)^ and as a result ?7™ can differ by 







, which 



makes the derived relation legitimate for odd m, but not for even m. We therefore prove that the fermionic SPT 
phases with are all equivalent to bosonic ones when m is odd. 



